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ABSTRACT 


The objective of this thesis was to develop a practical computer 
system for use of empirical data in the aerodynamic optimization of 
aircraft propellers. The system was designed for use with the IBM 
274lon-line computer terminal. This program provides instructions 
to the operator during execution, and allows interaction by the operator 
for input and alteration of data, and for program instructions. 

The Lockheed P-3C aircraft was chosen as the subject for test and 
evaluation of the program. The currently operational propeller of this 
aircraft was tested to compare the program's sredeeeton of aircraft 
performance against flight test information. An attempt was then made 
to select a propeller which would provide better performance under the 


same constraints as those imposed in design of the operational propeller. 
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I. INTRODUCTION 


In 1902, the Wright brothers, planning to adapt a theory used in the 
design of marine propellers, conducted a literature search in the Dayton, 
Ohio, public library. They found no theory, but only empirical data. A 
general theory of aircraft propeller performance was not forthcoming until 
1929 in a paper presented to the Royal Society of London by Goldstein, 
[Ref. 1]. This “vortex theory" developed by Goldstein has since been re- 
vised and amplified [Ref. 2], but is still cumbersome and somewhat ree 
accurate for use in real propeller design. 

Bang wind tunnel tests were conducted on aircraft propellers before, 
during, and immediately Poilieaing World War II, notably by Hartman and 
Biermann for the N.A.C.A. [Ref. 3]. Since that time the advent of the jet 
engine has caused a loss of interest in propeller research by all except 
-a few propeller manufacturers. 

Hence, even in 1970, aircraft propellers are still being designed by 
methods not much more sophisticated than those employed by the Wright 
brothers in 1902. A general theory exists which is limited in its appli- 
cability. Great volumes of empirical data, much of it useless in practice, 
has been generated. Still, however, the propeller designer basically 
searches through empirical charts for an aerodynamically optimum propeller 
configuration. 

Hamilton Standard Division of United Aircraft, probably the leader 
in the field of propeller design, has estimated sixteen man-hours of this 
type of work for each propeller design. Some use was being made of modern 
digital computers, however, to choose the best design after the range of 
possibilities had been narrowed by the manual process. The program. : sed 


were reportedly based upon Goldstein's vortex theory of 1929. 





The results of the extensive research and the general theories, in 
spite of the lack of sophistication, have been remarkable. The Wright 
brother's propeller, designed by a theory of their own, had an efficiency 
of about 66%. Modern propellers can be expected to exhibit efficiencies 
above 90%, in some cases. 

Judging from the history of the past seventy years, it seems reason- 
able to assume that empiricism will continue to dominate the propeller 
field in the foreseeable future. Therefore, it is the purpose of this 
thesis to make the procedure of designing an aerodynamically optimum pro- 
peller a little easier, and less time consuming, by employing the on-line 
digital computer systems now available. 

The result of this ae was a computer program, PROPOP, designed 
for use with the IBM 2741 on-line computer console. PROPOP asks ques- 
tions concerning the aircraft, engines, environment, constraints and opti- 
mization criteria. The required information is supplied by the Saeco 
through the communications console. 

PROPOP was designed with the intent that it contain all the instruc- 
tions needed by the operator to manipulate the program successfully. 

Parts II and III of this thesis serve as a user's manual, in case more 
detailed instructions on program use are required, or in the event the 
user wishes to modify the program in some way. Theoretically, no knowledge 
of FORTRAN or internal operation of the program is required for use of 
PROPOP, beyond the ability to use the IBM 2741 console (or similar system). 
Once execution of PROPOP begins, all instructions for its operation are 
provided by the program itself. 

The advantages of PROPOP lie in the integration of propeller efficency 


determination and aircraft performance evaluation. Present programs 
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apparently begin with assumed operating conditions for horsepower, alti- 
tude and velocity, then optimize propeller efficiency for these condi- 
tions. However, after this propeller is determined, it may be found 

that the aircraft does not actually operate at the assumed conditions. 
Consequently, the propeller is no longer optimum. The designer is then 
forced either to ignore this fact or to choose new assumed operating 
conditions and repeat the optimization process. In effect, PROPOP auto- 
matically chooses the latter of these two alternatives. Much flexibility 
is also incorporated into the system by allowing the operator to select 
the relative importance of the various aircraft performance characteris- 


tics, and by allowing him to impose constraints on aircraft performance. 
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II, PROGRAM ORGANIZATION AND OPERATION 


A. GENERAL ORGANIZATION 
As shown in Fig. 1, the complete program consists of three primary 


routines: 


ee Input Subroutine 
ae Optimization Subroutine 
Soin Performance Evaluation Subprogram 


Immediately upon execution of the program, propeller performance data 
are read through the off-line card reader. These data are internal to the 
program. A message is typed on the operator's console when all data have 
been read. A brief explanation of the program is then printed, followed 
by a list of steps the operator must perform. 

The Input Routine then begins. This routine asks questions of the 
operator about the aircraft, engines, etc., for which the propeller is 
to be optimized. The operator responds by typing in the required informa- 
tion and the program proceeds to the next question. 

PROPOP operates in two modes. If the operator completely specifies 
the propeller during operation of the Input Routine, the program operates 
in the evaluation mode. In this mode the program proceeds from the Input 
Routine to the Performance Evaluation Routine (dashed line, Fig. 1) and 
evaluates the aircraft performance with the given specified propeller. 
However, should the operator leave any or all the propeller parameters un- 
specified, the program operates in the optimization mode. In this mode, 
after completion of the Input Routine, the program proceeds to the Opti- 
mization Routine and optimizes for the unspecified parameters (solid line, 


Fig. 1). The Optimization Routine repeatedly calls upon the Performance 
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Evaluation Routine for information needed in the optimization procedure. 
At its completion, the Optimization Routine prints the optimum propeller 
parameters as well as the aircraft performance with this propeller. 

The next step in the procedure, regardless of the mode of operation, 
is evaluation of the output by the operator. At this point the operator 
may either terminate the program or he may return to the Input Routine 
and change any of the information that he has previously supplied. The 
usual action would probably be to change the constraints imposed on the 
aircraft or the propeller, or to change the priority given to any parti- 


cular portion of the aircraft's mission. 


Be INPUT ROUTINE 

The Input Routine, Fig. 2, consists of seven input steps: 

1. Aircraft data 

2. # Engine data 

3. Aircraft environment and anticipated performance 

A, Aircraft performance ‘constraints 

5. Propeller parameter constraints 

6, Optimization criteria 

7. End of routine (begin evaluation or optimization) 

Each of the first six steps requires the operator to furnish data to 
the program by typing it into his terminal. Step 7 is only a reference 
point within the routine so that the operator can refer to the end of 
the routine by calling for Step 7. The broken lines without arrows in 
Fig. 2 indicate that the operator may skip from one step to another, 
either forward or backward, within the Input Routine. For instance, 
Step 4 and Step 6 are not required if the program is to operate in the 


evaluation mode, so these steps may be omitted altogether. 
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FIGURE 2 
INPUT ROUTINE FLOW DIAGRAM 
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In each of the first six steps, PROPOP asks the operator for informa- 
tion required by the program. The operator then types these Pree into his 
console. At the end of each step, all the data supplied by the operator 
are printed on the terminal so that the operator may check it for ac- 
curacy. Should any data be incorrect, the operator simply returns to the 
beginning of that step and begins again to answer the questions asked by 
the program, typing in the correct data. 

The information called for by each step in the Input Routine is as 
follows: 

ie Step 1 (Aircraft Data) 

a. Design gross weight (lbs.) 

b. Wing area (sq. ft.) 

Cc. Number of engines 

d. Aspect ratio 

e. Parasitic drag coefficient (takeoff configuration) 
f. Parasitic drag coefficient (cruise configuration) 
g. Wing efficiency factor (e) 

h. Fuel capacity (lbs.) 

i. Takeoff speed (kts. TAS) 

Fae Step 2 (Engine Data) 

Engine data required by this step are determined by whether 
the aircraft uses reciprocating or turbine engines. 
a. Type of Engine (Reciprocating or Turbine) 
b. Reciprocating Engine 
(1) Takeoff brake horsepower 
(2) Engine RPM during takeoff 


(3) Brake horsepower during climb 
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(4) 
(5) 
(6) 
(7) 
(8) 
(9) 


Engine RPM during climb 

Brake horsepower during cruise 

Engine RPM during cruise 

Brake horsepower at maximum aircraft speed 
Engine RPM at maximum aircraft speed 


Fuel consumption (per engine) during cruise (lbs/hr) 


C. Turbine Engine 


(1) 
(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Engine operating RPM 

Shaft horsepower (per engine) at cruise altitude and 
cruise power, for anticipated cruise speed and for 

50 kts. above this speed 

Shaft horsepower (per engine) at maximum power setting 
(power for takeoff and maximum speed) and at 10,000 
ft. MSL 

Sea level shaft horsepower (per engine) for maximum 
power setting at zero, 100, and 400 kts. TAS 

Sea level jet thrust (per engine) for maximum power 
setting at zero, and 100 kts. TAS 

Sea level jet thrust for cruise power setting at zero 
and 100 kts. TAS 

Fuel ei on (lbs/hr) per engine for cruise power 
and altitude, at anticipated cruise speed and 50 kts. 
above this speed 


Limiting shaft horsepower for continuous operation 


All data requested for reciprocating engines should correspond to the 


referenced aircraft operating conditions. For example, "cruise horse- 


power” should be given for cruise RPM, cruise altitude and anticipated 


cruise speed. All conditions are specified above for turbine engines. 
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All engine data are for the installed engin’. Shaft horsepower 
(Item 4) for the turbine engine should be entered without regard to 
maximum allowable SHP due to gear box restrictions, etc. In other words, 
SHP's in Item (4) above may be greater than maximum allowable SHP in 
practice. This is necessary because of the interpolation process employed 
by the program. However, SHP's greater than the maximum allowable SHP are 
not actually used in calculating aircraft performance. 

3. Step 3 (Aircraft Environment and Anticipated Performance 

ae Takeoff altitude 

b. Cruise altitude 

Cc; Altitude for maximum speed calculations 

d. Anticipated speed for best climb rate (TAS) 

e. Anticipated cruise speed (kts. TAS) 

f. Anticipated maximum speed (kts. TAS) 

The altitudes specified above are those at which aircraft 
performance is calculated and are reference altitudes for much of the 
engine data requested in Step 2. The anticipated speeds are again refer- 
ence points for engine data, but are also starting points for the itera- 
tive process employed by the computer to calculate aircraft performance. 

4. Step 4 (Aircraft Performance Constraints) 

This step is required only for optimization of a propeller, not 
for simple evaluation of a given propeller. 

a. Maximum allowable takeoff distance (ft) 

b. Lowest allowatie maximum climb rate (ft/min) 

¢. Minimum allowable range at cruise conditions (n. mi.) 


d. Minimum allowable endurance at cruise conditions (hrs) 
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Ds 


e. Minimum allowable cruise speed (kts. TAS) 
f. Lowest allowable maximum speed (kts. TAS) 


Step 5 (Propeller Constraints) 


a. Maximum blade diameter for structural or ground clearance 
(ft) 

b. Number of blades (2, 3, or 4) 

c. Gear ratio (engine RPM/prop RPM) 

d. Blade diameter (ft) 

e. Activity factor (80-220 except for 2 blades, then 80-120) 

£. Integrated lift coefficient (.15-.7 or .3-.7 for 2 blades) 


Although this step is required for both the evaluation 


and the optimization mode of PROPOP, only the first item (maximum blade 


diameter} must be specified. Any or all of the other items may be left 


open to optimization by typing zero (0.) after the request for that parti- 


cular item. If any of these items is left unspecified in this manner, the 


program automatically operates in the optimization mode. Should all the 


items be specified, the program operates in the evaluation mode. The 


limits of the data currently available to the program are indicated above. 


6. 


Step 6 (Optimization Criteria) 


This step is required only for the optimization mode. 
a. Value of importance of short takeoff distance 

b. Value of importance of climb rate 

Cx Value of Pee cande of range 

ds Value of importance of endurance 

e. Value of importance of cruise speed 


f. Value of importance of maximum speed 


‘Values of Importance" are numbers from zero to ten indicating the 
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relative value of a particular aircraft performance characteristic to 

the overall aircraft mission. The absolute value of these numbers is 

not critical, only the relative values. Where more than one performance 
characteristic is given a non-zero value, the Optimization Routine stresses 
the characteristics in the order of their values, the highest value re- 
ceiving highest priority. 

The "Starting Point Routine" (Fig. 3) is designed for use in the 
optimization mode of PROPOP. Its purpose is to make an intelligent guess 
at the optimum propeller configuration before the actual optimization pro- 
cedure begins. This is necessary for three reasons. First, the Optimiza- 
tion Routine must begin somewhere and the closer it begins to the final 
optimum the better. Secondly the procedure must start within the range 
of propeller data available. The optimization procedure itself forces 
the propeller to lie within the — of data after the procedure begins. 
The third, and probably most important reason, is that the first propeller 
configuration chosen should be one with which the aircraft can fly. For 
example, it is difficult for the program to determine which propeller 
parameter to change for an improvement in cruise speed if the aircraft 
cannot even takeoff with the starting propeller. In such a case PROPOP 
attempts to optimize climb rate in order to find a propeller with which 
the aircraft will fly. This procedure is, however, subject to possible 
failure. 

The basis of the Starting Point Routine is the Aeroproducts Prelimin- 
ary Propeller Selection Chart, a generalized copy of which is shown in 
Fig. 4. The data in this chart are included in PROPOP and have the ef- 


fect of introducing an “experience factor" into the Starting Point Routine. 
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Explanation of terms in Fig. 3: 

BHP - brake/shaft horsepower; horsepower specified for maximum 
velocity for reciprocating engine, maximum continuous 
horsepower for turbine engine 

ALT - altitude at which maximum aircraft speed is to be 
evaluated 

RPM - engine RPM at maximum aircraft velocity 

B - number of propeller blades per engine 

AF - blade activity factor 

D - blade diameter 

CLI - blade integrated lift coefficient 

GR _ ake eae (engine RPM/prop RPM) 

EPA ~ effective propeller area (B x AF x p*) from Fig. 4. 

DMAX- maximum allowable propeller diameter 


J  - advance ratio (v/nd) 


TAS, ft/sec 


< 
" 


n = propeller RPS 


oS 
if 


Propeller diameter, ft. 


C. OPTIMIZATION ROUTINE 

The optimization procedure is based on a pattern search technique 
developed by Hooke and Jeeves [Ref. 4, 5]. The subroutine DIRECT, in- 
corporated into the IBM System/360 Subroutine Library in 1966 was used 
as the basic optimizer. It is this procedure which is outlined in Figs. 
5 and 6. This subroutine is, however, an unconstrained minimization pro- 
gram. For this reason the routine had to be heavily modified so that it 


could handle the highly constrained propeller optimization problem. 
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Basically, pattern search is a numerical trial-and-error procedure 
closely akin to the better known "gradient search" methods (Ref. 6]. A 
fundamental difference, however, is that pattern seers. unlike gradient 
methods, does not attempt to follow the path of steepest descent to the 
minimum of the function value. It is this very property which often foils 
gradient methods when they encounter a "ridge" or "steep valley" in the 
objective function. The pattern search begins by exploring the region of 
the starting (base) point by incrementing each variable separately in 
search of a function value which is better (lower in this case) than the 
base point (Fig. 6). After having explored in all the separate variable 
directions, the program will have normally found a Bega: function value 
than the base point value. It then assumes that further decreases in the 
function value may be made in this same direction from the base point 
(Fig. 5). The program has now established a "pattern" and will continue 
to make pattern moves in this direction until a move fails to produce an 
improvement. At this point another exploratory eee is performed and 
the procedure begins again. If all exploratory steps also fail to pro- 
duce an improvement, the step size is decreased to increase the defini- 
tion of the functional surface. The program terminates when the step 
size falls below a prescribed value. 

Modification of this subroutine included built-in “stops' to insure 
that the search procedure does not leave the bounds of the available data. 
Other modifications allowed the operator to specify any of the variables 
so that the search procedure does not change these variables. The speci- 
fied variables are known as equality constraints because they are always 


equal to a set value. Constraints which specify a range or a limit of 
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acceptable values of the variables, such as the data limits, are known 
as inequality constraints. 

Another set of important inequality constraints are the aircraft 
performance constraints. These constraints were accomplished by modi- 
fications of the Performance Evaluation Routine rather than the opti- 
mization procedure itself. if at any time the propeller variables being 
tested by the Optimization Routine result in aircraft performance which 
is unacceptable by the operator's prescribed standards, a false function 
value is assigned for that trial. This false function value is designed 
so that it decreases very rapidly as the "boundary of acceptance” of the 
aircraft performance is approached. 

In order to understand the operation of this technique more clearly, 
it is necessary to have an explanation of the way in which a functional 
value of aircraft performance is derived. The technique employed by 
PROPOP is as follows. For a given propeller configuration the nO REOrTence 
of the aircraft is evaluated by the Evaluation Routine. A functional 
value (which has been derived by the author) is then assigned. This 


value has the form: 


/ 
Fos 2, {(vOr) Pee Nee] 

where 

F is the functional value 

i is an index denoting one of the performance criteria (takeoff 

distance, maximum speed, etc.) 
VOT) | is the value of importance assigned by the operator to the ee 
performance criterion 

P. is a normalized measure of the ie performance criterion which 


i 


th : 
always decreases as the i performance criterion improves 
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(P. is of order 1) 

N is simply an indicator which equals zero if the ‘a criterion 
is acceptable and equals one if the ike criterion is unaccept- 
able 

P is another measure of the Aes performance criterion which de- 
creases as the has criterion improves, but P is of order 100 
or greater. 

The Optimization Routine secks to reduce the value of F, so as long as 
each performance criterion meets the constraints imposed by the operator, 
the routine attempts to decrease the value of P. and consequently im- 
prove the ie performance criterion. Because all P. are of order one, 

a high value of (VOT) . gives more weight to the ae criterion in the 

overall function. However, if cone or more of the constraints is not 


V6 
satisfied, then N, equals one, which allows P. to dominate the other terms 


i 
of the function. The optimizer immediately goes to work on improving the 
gy criterion until the boundary of acceptance is crossed, at which time 


the routine continues to optimize the P, values giving preference to 


i 
higher values of importance. 
Five variables are recognized by the Optimization Routine. These 

are: 

1. Number of blades 

2. Gear ratio (engine RPM/prop RPM) 

cP Propeller diameter 

4. Blade activity factor 

Sie Blade integrated lift coefficient 


Of these five, the first is a discrete variable, while the other four are 
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continuous variables. In other words, there can be only an integer 
number of blades, but the other variables may take on fractional values. 

Discrete variables present a problem for numerical optimization 
routines. There are two widely used methods of treating this problem. 
Probably the most obvious method is to search for a local optimum with 
each possible value of the discrete variable, separately, finally choosing 
that value for which the best local optimum is found. The second is to 
treat the discrete variable as a continuous variable, allowing it to 
have fractional values, and then rounding off to the nearest integer 
value. 

Both methods have advantages and disadvantages. The first approach 
would obviously be very time consuming if the discrete variable had many 
possible values. However, it would eventually find the correct optimum. 
The second method would be faster, but the rounding off procedure could 
easily be fooled. In practice, with a large range of possible erie’ 
for the discrete variable, a combination of the two methods may be used. 
The continuous variable approach could be employed to narrow the choice 
to a few values of the variable, each of which could then be searched 
separately. 

In the propeller optimization problem, data were available only for 
2, 3, and 4-blade propellers, so there are only three possible values of 
the discrete variable. The individual search approach was therefore 
employed in the Optimization Routine. Figure 7 outlines this procedure. 
Unless the number of blades is specified by the operator, the routine 
first conducts a search with a 4-blade propeller, followed by a 3-blade 


search. If the 3-blade optimum is worse than the 4-blade optimum, the 
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search procedure ends, assuming that a 2-blade propeller would result 
in a still more inferior performance. However, should the 3-blade pro- 


peller outperform the 4-blade prop, then a 2-blade search is conducted. 


D. |§ PERFORMANCE EVALUATION ROUTINE 

The purpose of the Performance Evaluation Routine is to calculate 
aircraft performance for a given propeller configuration. This is ac- 
complished by use of propeller performance maps found in Ref. 7. These 
charts give static thrust data and efficiency data for propellers of 
varying number of blades, activity factor, and integrated lift coeffi- 
cient. These charts were reduced to data points and are fed into the 
program on punched cards. 

According to Hamilton Standard Division of United Aircraft, the 
authors of Ref. 7, these data were derived semi-empirically: partially 
by wind tunnel data and partially by computer computation using a modi- 
fied vortex theory originally developed by Goldstein in 1929. The 
technique employed by the PROPOP routine is to repreduce these perform- 
ance charts by Lagrangian non-linear interpolation of second order. It 
was found that when each efficiency map was reduced to 42 data points, 
the accuracy with which this technique could reproduce the full chart 
was normally within 1%. This was considered quite adequate, as the data 
are assumed only to be accurate to within 1 - 2%. 

Another approach to this calculation might have Beek to start with 
vortex theory and let the computer compute the efficiencies just as 
Hamilton Standard has done. In this way the program would not have been 
limited to the range of data presented in the charts. However, since 


data were already available covering a rather substantial range, the 
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added complexity of the program and the additional computer time which 
would be required to reczlculate this information were considered too 
great a burden to endure for the amount of aesthetic satisfaction which 
might have been obtained. 

Corrections for compressibility and nacelle blocking effects are 
presented in graphical form in supplements to Ref. 7. These charts were 
reduced to data points in the same manner as the efficiency charts and 
were incorporated into the PROPOP data bank. The accuracy of these cor- 
rections is somewhat less than the accuracy of the original data, but 
they are useful in determining trends. The general operation of the 
Evaluation Routine is outlined in Fig. 8. 

Turbine engines Present a problem when it is necessary to calculate 
their performance under a wide range of operating conditions. In a pro- 
gram such as this, it was considered important to reduce the amount of 
data required of the operator to a minimum. It was therefore necessary 
to make assumptions on the variation of engine performance with power 
setting, altitude, and velocity. The assumptions made are as follows: 

1 Shaft horsepower at a given power setting is assumed to 
vary linearly with altitude and non-linearly (second-order) with velocity. 
2. Jet thrust at a given power setting is assumed to be con- 
stant with altitude, and to vary linearly with velocity. 
3 Fuel consumption at cruise power and altitude is Pere 
to vary linearly with velocity. 
| 4. Constant engine RPM and constant gear ratio. 
5. emeeee. maximum speed, and climb rate are calculated at 
same power setting ("maximum power"). 
6. Range, endurance, and cruise speed are calculated at 
"eruise power" setting and cruise conditions. | 
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With the above assumptions and the data furnished to the program by 
the operator, engine performance is calculated at all necessary operating 
conditions. 

L. Calculation of Takeoff Distance 

a. Assumptions 

a C= C,+ C, 


where: C, = aircraft parasitic drag coefficient in takeoff configuration 
Pp 


Cy = induced drag coefficient 
i 
2 
EC = 1C /t AR e 
ws, L 
where: AR = wing aspect ratio 


@ = wing planform efficiency factor 


C. = aircraft lift coefficient ° 
, 2 
= ai _ a= + 
cy trust = Teast [Lo Tsrarlye= Gras? BV 
where ae static thrust (zero velocity) 
~. = thrust at liftoff | 
\J/ = instantaneous velocity (TAS) 


Vo = liftoff velocity (TAS) 


. IT. O re VA Mo 


Thrust is assumed to vary parabolically with velocity during takeoff 


OO 


run. 
b. Calculations 
acceleration = a = F/M = summation of forces/aircraft mass 


F = thrust - friction drag - aerodynamic drag 


van 0 QW LEAVES [6p G_ 


where: b- = coefficient of friction between runway and wheels 
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W = takeoff weight of aircraft 
ia = aircraft lift 


Pp = air density 








S = aircraft wing area 
4 
OL: = { 2 | 2 
— —_— 4+ — 
D=p[W-sPV*SC |+zpV S|¢,+ Tare | 
In order to find the optimum C, for minimum drag: 
ap VPS oes 
dC PP Tare toot 
= + wllare 
tae im 


Aircraft is assumed to maintain a constant C. (Cc, ) during takeoff run. 
opt 


a- S17, t8 v= p|w- +pv?s Cah tpv' S c} 
os a 1 {(Torar Wt [Bt 3p SKC, ~ Cyl] V4} 
Takeoff distance (x, >): 


va a/v gyav= fy 


ot fay 
A+z\v* 





where: A = Tstat pW 
“(Cho Wsrat Ete 2 TPs S (nC, Cy) 


: air density at sea level 
T =P/ Ps. 
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Further assumptions: 

OD) UL 025 

(2) no wind 

(3) no runway slope 

(4) standard temperature 

(5) Oo = exp (-altitude/28900) 

This procedure is explained in greater detail in Ref. 8. 
Takeoff distance can vary greatly with pilot technique, temperature, 
etc., but this is of little concern to the optimization of a propeller. 
In general, a good takeoff propeller under the assumed conditions would 
also be a relatively good takeoff propeller under other conditions. 
Relative performance of one propeller configuration to another is the 
important factor. 

Ze Calculation of Endurance and Range 

a. Endurance 
Endurance is calculated at cruise speed, cruise 
power and altitude on the basis of a 45 minute reserve. 
Endurance = (fuel capacity/fuel consumption)-45 minutes. 
Fuel consumption is assumed constant. 
b. Range 
Range is also calculated at cruise speed, 


cruise power and altitude on the basis of a 45 minute reserve. 


Range = Endurance x V 
ave 


36 





where: Neen = average velocity during flight or 
V +2 AV 
AV = cruise velocity with empty fuel tanks minus crulse 
velocity with full fuel Ces 


AV is estimated as follows: 





1p =liV= DV _ y = Pavan .550 
Avail 550 # £550 D 
where es 
_ aed 2 
Pont Gale zpV Ss 
and —_ 2W 
C, pVrs 
; S5018 i a WE 
VS LOo+t+—Se 
_T O pv" Stare 


(0° S*TareC, Vv" IOOFP p Siar eV +4W* = O 


AVAIL 
Taking the derivative with respect to\/ : 
w= alma 
p Stheele75Fp — p SC, V° 


AVAI 





or approximately 





2W AW / At 
= ——= — = | At 
AV Reeaa it Ren qo Cy, “A 
where: 
AV = — ay 


so initial 


At = endurance 


—AW/At = fuel consumption rate 


\V = cruise velocity (V_) 
cr 
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W = W, GROSS . o> Weyer = average weight 
Accounting for increased velocity due to loss of fuel weight typically 
increases range by 2-3%. 

Range and endurance are calculated at cruise speed and conditions 
for seweral reasons. If given conditions were not specified, these 
performance characteristics would have to be determined by an iterative 
process of varying the engine power setting to determine the best setting 
and the corresponding aircraft velocity for best endurance. Such a pro- 
cess, although conceptually straight forward, would require more data 
than are presently Paruited by the Input Routine. For instance, engine 
fuel consumption would have to be known throughout the entire range of 
engine operation. The iteration process would also require extensive 
computational time. In fact, such a procedure would probably double the 
time required by the entire Evaluation Routine. These penalties were 
eonsidered too great to justify the computation. 

3. Calculation of Maximum Velocity and Cruise Velocity 
Calculation of maximum aircraft velocity and cruise veloc- 
ity are essentially identical and the computational procedure is outlined 
in Fig. 9. Basically the procedure involves locating the intersection of 
the power available curve and the power required curve (both vs. velocity) 
by iterating velocity [Ref. 8]. 
The symbols in Fig. 9 are as follows: 
A + step size or velocity increment 
N1,N2 - indicators which serve as cues for the routine to 
branch in the proper direction at certain points 
rl - thrust horsepower required for level flight at given 


conditions 
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velocity step size was set at .5%.of the anticipated climb velocity, so 
this is the accuracy of the resulting climb velocity. However, at the 
maximum climb speed, excess thrust horsepower normally varies very slowly, 
so the accuracy of the actual rate of alae calculation should be much 
better than .5% given perfect information. 
The symbols used in Fig.10 are as follows: 
~ step size (.005 times anticipated climb speed) 
N| - indicator which serves as cue for a switching point in 
the routine 
VV - velocity of aircraft used in calculating excess thrust 
horsepower (TAS) 
AHP ,AHP - excess thrust horsepower at two consecutive velocity 
trials 
HP - required thrust noe: for level flight under given 
conditions 
HP - available thrust horsepower at given conditions 
Woross - aircraft gross weight 


R.C. - rate of climb (fpm) 
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III. PROGRAM USAGE 


A. OPTIMIZATION OF PROPELLERS FOR GIVEN AIRCRAFT 

The first step in optimization of a propeller by use of PROPOP is 
the gathering of the required data, as listed in Section II B, on the 
aircraft, engines, and operational requirements. It is important that 
all engine horsepower (and jet thrust, if turbine powered) information 
used is installed engine data. Specifications as to operating conditions 
of the engine must be followed exactly or erroneous calculations will re- 
sult. 

One of the features of PROYOP is that the horsepower limit specified 
by the operator is automatically imposed on turbine engines. It is there- 
fore unnecessary (and undesirable) for the operator to impose this limit 
on the horsepower data typed into the program (see Section II B.2.c.). 
PROPOP recreates engine power curves by Bipeiee given data points with 
either second order curves or straight lines. These curves could be 
drastically altered if the operator imposed such artificial constraints 
on the data. It may be necessary, therefore, to extrapolate engine power 
curves into regions where the engines do not actually operate in order 
to provide the program with smoothly varying curves by which trends may 
be calculated. 

Steps 1, 2, and 3 of the Input Routine are required for operation of 
PROPOP in the optimization mode. Step 4 (Aircraft Performance Constraints) 
may be omitted altogether. In such a case no constraints are imposed on 
the ee performence. Should the operator desire to specify some, 
but not all the performance constraints, he may indicate those to be left 


unspecified by assigning the value zero (0.) to those constraints. 
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Step 5 (Propeller Constraints) is required, however any constraint, 
except "maximum blade diameter," may be left unspecified by assigning a 
value of zero (0.) to that constraint. The range of data acceptable for 
this step is given in Section II B.5. 

Step 6 (Optimization Criteria) is required in the optimization mode. 
In this step the operator specifies the aircraft performance character- 
istics for which he requires an optimum propeller. He may specify one or 
more of these characteristics (Section II B.6) by assigning "values of 
importance" from zero to ten. The highest values of importance receive 
the highest priority in the Optimization Routine. It has been found that 
the best technique for assigning these values of importance is to specify 
only one, giving zero (0.) as the value of the others. An attempt has 
been made to insure that the resulting aircraft performance closely re- 
flects the relative values of importance assigned by the operator (in the 
case of several specified values of importance), however this is not always 
the case. The operator should specify only one characteristic, anck as 
cruise performance, for the first trial. If the resulting propeller con- 
figuration gives other performance characteristics which the operator 
considers undesirable (such as large takeoff distance), then he may re- 
turn to the Input Routine (Step 4) and specify a maximum allowable take- 
off distance or other constraint. 

After much testing of PROPOP it was found that most of the objective 
functions which must be optimized exhibit a "broad bottom" phenomenon. 
This means that in the region of the desired minimum fairly large changes 
in the variables may result in very small changes in the function being 
optimized. As a result of this situation there are many combinations of 


propeller configuration variables which yield nearly the same value for 
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propeller configuration variables which yield nearly the same value for 
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a particular aircraft performance characteristic. However the perfor- 
mance characteristics not considered in the optimization may vary greatly 
with each configuration. For example, two different combinations of gear 
ratio, diameter, and activity factor may yield the same optimum cruise 
speed, but takeoff distance for one combination of these variables may 

be twice that for the other configuration. 

Due to "noise'' in the data, this broad bottom phenomenon may cause 
problems of apparent multimodality, or multiple optima. The operator may 
be surprised and somewhat disillusioned if he optimizes twice for the same 
performance characteristic and PROPOP gives two very different propeller 
configurations. However it should make no difference which configuration 
the operator chooses, sitrice they both provide the same performance. This 
problem may be greatly reduced by specifying constraints on several per- 
formance characteristics, as mentioned above, while optimizing for a 
single characteristic. 

The PROPOP output permits the operator to follow the program's pro- 
gress through the optimization procedure. After completing the search 
for a particular number of blades, PROPOP gives the number of function 
evaluations performed. This is the number of trial solutions attempted 
in the search. Five hundred (500) solutions are the maximum permitted. 
If 500 function evaluations are performed, PROPOP gives zero (0.) for 
the value of this number. This is aii indication to the operator that 
an optimum solution has not been reached, and if more accuracy is desired 
the limiting number of evaluations must be raised above 500. This may be 
done simply by changing the next to last calling argument in the "CALL 


DIRECT" statement, which is statement number 15, near the end of the 
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INPUT SUBROUTINE (Appendix A). However, such a problem is highly un- 
likely, as the use of more than 200 evaluations is a very rare occur- 
rence. More information on the optimization subroutine (DIRECT) may 
be found in the comments preceeding that subroutine (Appendix A). 

For the direct search optimization procedure (Section II C.), best 
results were found by starting with a rather large step size and de- 
creasing it slowly until the desired accuracy was obtained. Consequently, 
the beginning normalized step size (DELCAP) was set at .816. The step 
length is decreased to one fourth of its current value (RHO = .25) when- 
ever finer resolution is required. The minimum step size used is .012/. 
This value has proven to be sufficiently small, and further reduction 
rarely changes the results of the search. However, should the operator 
so desire, any of the above values may be changed by consulting the 
comment section of SUBROUTINE DIRECT and then altering the calling argu- 
ments in the "CALL DIRECT" statement previously mentioned. DELCAP is 
changed by altering statement 8 in SUBROUTINE INPUT. The current values 
of these arguments are printed in the PROPOP output before each search. 

DIRECT operates by assigning the same step size to all variables. 
This is obviously a disadvantage when the variables have vastly differ- 
ent ranges, such as zero to infinity, and .15 to .7. Consequently, it 
was arranged that DIRECT would deal only with normalized variables rather 
than the actual ones. The actual variable values are therefore divided 
by an assigned step size before being used by DIRECT. The step sizes 
assigned are as follows: 

diameter.........-1 ft. (actual value of variable) 


PGdmeratlO.....ce..+. «2 (CSTG) 
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activity factor ... 90 (CSTA) 


integrated C we, (CSTL) 


Lott 
So the beginning step size for activity factor, for example, is 73.4 

(or .816 x 90). In selecting these step sizes, an attempt was made to 
have the contours of the objective function be as nearly circular as 
possible near the optimum. In other words a change of one step size in 
any of the variables should have approximately the same effect on the 
propeller performance. Such a condition makes optimization more effi- 
cient. These step sizes, except that for diameter, may be altered by 

the operator by changing the values of CSTG, CSTA, or CSTI which are 
found near the beginning of PROPOP (Appendix A). 

Much effort has been placed on restraining the optimization search 
to the limits of the data available to PROPOP. This is a relatively 
simple task for variables such as activity factor and integrated lift 
coefficient because of the well defined limits to these variables in the 
data. It is very much more difficult, however, for power coefficient 
(c)) and advance ratio (J). The difficulties arise for several reasons. 
Power coefficient, for example, has no well defined limits in the data; 
the limits being different for each chart in Ref. 7. Another problem 
is that the search procedure does not actually fix c, or J. Instead it 
sets other propeller parameters and later calculates oe and J while 
determining aircraft performance. By that time errors may already 
have occurred. Because of these complexities, no attempt was made to 
constrair. oe Advance ratio, J, is constrained rather crudely to an 
upper limit of three (3). The value of J is computed (V/ND) assuming a 
velocity which is the anticipated maximum velocity specified by the oper- 


ator. If this value is greater than three, PROPOP alters either gear 
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ratio, to change propeller RPM (N), or diameter (D), so that J equals 
three. If the actual maximum velocity is greater than that assumed by 
the operator, however, J will still lie outside the available data. A 
small overshoot in oa or J is usually not fatal, as the extrapolation 
employed by PROPOP is usually fairly accurate close to the data limits. 
Large overshoots, however, can be disasterous, yielding highly inaccurate 
results (possibly efficiencies less than zero or greater than 100%). 

If the value of J corresponding to maximum speed should be very 
close to three after an optimization, there is a good possibility that 
the propeller diameter or gear ratio have been altered during the pro- 
cedure to restrain the search to the data limits. In this situation it 
is sometimes possible to obtain a better optimum by increasing the assumed 
maximum speed (and changing the appropriate corresponding engine data) 
aC the risk of overshooting the bounds of the data. 

Due to the process employed in imposing limits on propeller diameter, 
L? PROPOP occasionally gives an optimum 


propeller with one or more of the parameters very close (well within one 


activity factor and integrated C 


etal step size) to its limiting value. In such a case, a slightly 
BeGeer optimum may sometimes be found by constraining that parameter to 
its limiting value. The operator should not be too surprised, however, 
if the results of this maneuver yield a propeller with the other variables 
radically changed from those of the first “optimum" propeller due to the 
multimodal phenomenon discussed previously in this section. In any case 
the aircraft performance should not vary more than one percent. 

Output provided by PROPOP includes the following information on air- 


craft performance with the optimized propeller: 
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Takeoff 


Takeoff distance 

Static Power Coefficient 

Static Thrust Coefficient/Power Coefficient (C,/C)) 
Propeller efficiency at liftoff 

Blade angle (.75R) at zero velocity 

Power coefficient at liftoff (CP) 

Advance ratio at liftoff (J) 

Effective Mach number at liftoff (M(EFF)) 


Compressibility correction applied to prop efficiency 
(FT) 


Blocking correction applied to prop efficiency (FB) 
Brake/shaft horsepower at liftoff (BHP) 


Jet thrust (turbine) (TJET) 


Maximum Speed 


a. Maximum speed 

b. Blade angle (.75R) 

c. Prop efficiency 

d. Power coefficient (CP) 

e. Mivence ratio (J) 

£. Effective Mach number (M(EFF)) 

g. Compressibility correction applied to prop efficiency 
(FT) 

h. Blocking correction applied to prop efficiency (FB) 

i. Brake/shaft horsepower (BHP) 

j. Jet thrust (turbine) (TJET) 

Cruise 

a. Cruise speed 

b. Same values as b.-j. for maximum speed 
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Cz Range 
a, Endurance 
4. Climb 
a. Maximum rate of climb 
Be Climb speed 
OF Same values as b.-j. for maximum speed 
Effective Mach number is the aircraft Mach number corrected by 


a factor involving blade integrated C Effective Mach number, compressi- 


L 
bility corrections (FT), and blocking corrections (FB) are explained in 
supplements to Ref. 7. Listed values of prop efficiency are for correc- 
tions already applied. The blocking corrections are ie percent below 
those given in supplements to Ref. 7 to account for interference and 
slipstream effects not included in the propeller efficiency data. This 
additional factor was suggested by Lockheed California Company [Ref. 9]. 
For takeoff performance, where ground effect is present, this additional 


blocking factor is omitted and only that factor given in supplements to 


Ref. 7 is applied. 


B.  EWALUATION OF A GIVEN PROPELLER CONFIGURATION 

As explained in Section II B., PROPOP may be operated in the evalua- 
tion mode by specifying all blade parameters in Step 5 of the Input 
Routine. Again, as in the optimization mode, Steps 1, 2, 3, and 5 are 
required. However, Steps 4 and 6 are omitted. 

Output data on aircraft performance is identical with that des- 
cribed in part A of this section for operation in the optimization mode, 
except of course there is no output corresponding to the optimization 


procedure itself. 


50 





IV. EXAMPLE PROBLEM: LOCKHEED P-3C 


The Lockheed P-3C aircraft was chosen as the test case for the 
PROPOP system. Choice of this aircraft was based on its use by the Navy 


and availability of data on the aircraft and propellers. 


A, EVALUATION 

As nearly as could be determined from information provided by the 
propeller and aircraft manufacturers, the P-3 propeller was optimized for 
a low power cruise (about 60% NRP) at an altitude of 15,000 feet, and at 
a speed of 325 kts. (TAS). It appears, however, that the aircraft will 
not fly at these conditions at full gross weight (135,000 lbs.). There- 
fore, a reduced weight of 105,000 pounds was used in calculations of 
cruise performance. Even at its empty weight of 66,211 pounds, the air- 
craft will not attain 325 kts. It therefore seems possible that the 
propellers were designed for conditions at which the aircraft cannot 
operate. 

Data used in the P-3 performance evaluation is given in Table I. 
PROPOP calculates takeoff distance, maximum speed, and rate of climb at 
the same power setting (‘maximum power"). It was therefore necessary to 
conduct two evaluations. Takeoff distance was calculated using takeoff 
power (1970°F T.1I.T.), while maximum speed and climb rate were calculated 
at military power (1920°F T.I.T.). A third evaluation was necessary for 
cruise speed, range and endurance, inasmuch as a reduced gross weight was 
employed. | 

It may be noted in Table I that several of the shaft horsepower 
figures are considerably above the 4500 SHP gearbox limit. This necessity 


wag discussed in Sections II B.2. and III A. 
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TABLE I 


P~-3C DATA SHEET 


AIRCRAFT DATA: 
Gross Weight 


Wing Area 

Number of Engines 

Aspect Ratio 

Parasitic Drag Coeff., C (takeoff) 

$9 st 08 D e 
(cruise) 

Wing Efficiency Factor (e) 

Fuel Capacity 


Liftoff Speed 


ENGINE DATA: 


Engine RPM 

SHP (assumed cruise power, 
altitude, and speed) 

SHP (assumed cruise power, 
altitude, and speed + 50 kts.) 
SHP (max. power, 10,000 ft. MSL, 
assumed max. speed) 

S.L. SHP (max. pwr., 0 kts.) 


Seieone (" Ee ,100 kts.) 
S.L. SHP (" "400 kts.) 
S.L. Jet Thrust (max. pwr., 0 kts.) 
19 08 18 ay 18 , 100 kts.) 
ae 4 (cruise " , O kts.) 
8 iT) 8 Cc" 7] , 100 kts.) 


Fuel Flow (assumed cruise power, 
altitude, and speed) 

Fuel Flow (assumed cruise power, 
altitude, and speed + 50 kts.) 


Limiting SHP 


ENVIRONMENT & ANTICIPATED SPEEDS: 


Takeoff Altitude 

Climb ‘Altitude 

Cruise Altitude 

Max. Speed Altitude (assumed) 
Anticipated Climb Speed 

" Cruise " 


oe Maximum 
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135,000 (lbs.) 
(105,000 for cruise) 
1300 (sq. ft.) 

4 


Use 
-05 (10 AERO 
-0247 15D PYLONS) 
Ss: 


59,370 (lbs) 
(29,370 for cruise) 
136 @ (kts. TAS) 
121 @ 105,000 lbs. 


13,820 
1755 1450°F 
geal fad 

1861 @ 25°C 

4453 (T.0. @ 4637) 

4305 (T.0. @ 4527) 

4486 (T.O. @ 4712) 

5560 (T.O. @ 5790) 
0°C 

627 (T.O. @ 639) lbs. 

520 (T.0. @ 533) lbs. 

462 lbs. 

353 lbs. 
1450°F 
Teil. 
@ 25°¢ 


1035 1lbs/hr 
1065 lbs/hr 


4500 (gearbox limit) 





Results of this performance evaluation and comparison with avail- 
able flight test data are presented in Table II. Good correlation was 
obtained between the predicted and the flight test values. The largest 
error was obtained in takeoff distance (about 4%). The predieced value 
may be slightly conservative, but considering the many variables con- 
trolled by the pilot during takeoff, this is not thought to be a signi- 


ficant error. 


TABLE II 
COMPARISON OF FLIGHT TEST DATA AND PREDICTED 


PERFORMANCE OF P-3C AIRCRAFT 


PERFORMANCE FLIGHT TEST PREDICTION 
Takeoff Distance 4300 ft. 4459 ft. 
Maximum Speed 386 kts. 388 kts. 
Cruise Speed * 287 kts. 
Range * 1742 n. mi. 
Endurance * 5.9 hrs. 
Rate of Climb 2010 fpm. 2086 fpm. 


*Flight test data not available for chosen conditions. 


B. OPTIMIZATION OF P-3 PROPELLER 

As an example, the P-3 propeller was optimized for the low power 
cruise discussed in part A of this section. It was assumed that the 
gear ratio (13.54:1) was determined by the engine manufacturer and was 
therefore constrained. The maximum allowable propeller diameter for 


structural clearance of the P-3C is 13.5 feet. Takeoff distance was 


constrained to under 2500 feet at a gross weight of 105,000 pounds. 
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A sample computer printout of the input sequence and the PROPOP out- 
put is found in Appendix B. The optimum propeller characteristics may 
be compared with the P-3 propeller {Ref. 10] in Table III. They vary 
only slightly, however the performance of the optimum propeller exceeds 
that predicted for the operational P-3 propeller in every respect. The 


performance comparison is given in Table IV. 


TABLE III 
COMPARISON OF OPERATIONAL P-3C PROPELLER 


AND SAMPLE OPTIMUM PROPELLER 


P-3 OPTIMUM 
Number of Blades 4 & 
Gear Ratio 13.54:1 13.54:1* 
Diameter «13,5 fe. 13.5 £t** 
Activity Factor 163 220%**x 
Integrated C . 286 as a 


L 


*Assumed given 
*«kMaximum diameter for structural clearance 
*kkUpper limit of data available 
TABLE IV 
PREDICTED PERFORMANCE COMPARISON BETWEEN 


OPERATIONAL P-3 PROPELLER AND SAMPLE OPTIMUM 


PERFORMANCE P-3 OPTIMUM 
Takeoff 4459 ft. 3875 ft. 
Maximum speed 388 kts. 389 kts. 
Cruise speed 287 kts. 290 kts. 
Range 1742 n. mi. 1853 n. mi.* 
Endurance 2.9 hrs. | 6.2 hrs.* 
Climb 2086 fpm. 2202 fpm. 


*Based on 29,370 pounds of fuel. (105,000 pounds gross weight) 
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V. EVALUATION OF PROJECT SUCCESS 


It is felt that the original objectives as previously stated have 
been realized successfully. A practical, relatively easy-to-use com- 
puter program has been developed for use of empirical data in the pro- 
cess of optimum propeller selection. PROPOP has been shown to predict 
aircraft performance with acceptable accuracy and to improve on existing 
hardware design. 

Time required from collected data to final propeller selection has 
been reduced from the sixteen man-hours necessary with a manual system 
to one half man-hours using PROPOP. A 97% reduction in required time 
with acceptable accuracy is considered sufficient incentive for further 


investigation into the feasibility of this system. 
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VI. PROBLEMS ENCOUNTERED AND AREAS FOR FUTURE RESEARCH 


Several problems were encountered in the development of the PROPOP 
system which, although apparently solved, may cause future difficulties. 
One such area is the constraint procedure employed during optimization, 
especially concerning Oe and J as discussed in Section III A. 

Further improvements could be made in the prediction of turbine 
engine performance. It appears possible that excessive data are presently 
required and the resulting predictions are not as accurate as might be 
hoped. Some simple functional relations for performance versus throttle 
setting, altitude, and velocity would be highly desirable. 

More flexibility and a more sophisticated calculation of range and 
endurance would also be desirable. At present these characteristics are 
calculated rather crudely at cruise conditions. 

The optimization procedure could also bear much more work. A 
system by which final aircraft performance relates more closely to assigned 
“values of importance" would be a significant improvement. The pattern- 
search algorithm upon which the optimization depends should be more 
thoroughly investigated to determine if another algorithm might be prefer- 
able for the propeller optimization problem. 

The "broad bottom" phenomenon observed, which often leads to multi- 
ple optima, should be investigated further. The possibility exists of 
rescaling the variables in an attempt to alleviate this problem. Other 
approaches might also be considered. 

The determination of the optimum relationship between the propeller 


variables for producing near-circular contours near the optimum is an 
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area which might be considered. The optimum initial step length and the 
best rate of reducing this length could also lead to improvements in 
program efficiency. 

If accurate theoretical methods of predicting propeller performance 
could be substituted for the empirical data and interpolation procedure 
employed by PROPOP, there exists the possibility of much improvement. 
Data could be eliminated, accuracy might be improved, and available data 
limitations could be removed which would eliminate many of the constraint 
troubles. Such a substitution might even solve the multiple-optima 
problem. 

PROPOP in its present form is peaiily only a beginning, with vast 
room for improvement. For instance, there seems to be no reason why a 
similar method could not be developed for optimal design of an entire 


aircraft. 
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APPENDIX A 


COMPUTER PROGRAM 
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APPENDIX B 


SAMPLE COMPUTER OUTPUT 
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